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Abstract. Vertex models with quantum group symmetry give rise to 
integrable cellular automata at g = 0. We study a prototype exam- 
ple known as the periodic box-ball system. The initial value problem 
is solved in terms of an ultradiscrete analogue of the Riemann theta 
function whose period matrix originates in the Bethe ansatz at g = 0. 

1. Introduction 

The periodic box-ball system ^1 ^] is a completely integrable one-dimensional 
cellular automaton. Its dynamics is described as a motion of balls hopping exclu- 
sively along the periodical array of boxes having capacity 1. The system is identified 
with a solvable vertex model 2 associated with quantum affine algebra Uq^slz) at 
(7 = 0, where the fusion transfer matrices Ti,T2,. . . yield a commuting family of 
deterministic time evolutions. 

In IIQI , the initial value problem of the periodic box-ball system is solved by an 
inverse scattering method. It is done by synthesizing the combinatorial versions 
of the Bethe ansatz at g = 1 [S] and q = The action-angle variables are 
introduced by generalizing the rigged configurations (g = 1) up to some equivalence 
specified by the string center equation {q = 0). It enables one to determine the 
time evolution T^p) of any state p by an explicit algorithm whose computational 
steps are independent of the time t. 

The Bethe ansatz approach ^U] captures several characteristic features in the 
quasi-periodic solutions of soliton equations For instance, the original nonlin- 

ear dynamics becomes a straight motion of the Bethe roots (angle variable) which 
live in an ultradiscrete analogue H2.6|l of the Jacobi variety. 

In this paper we exploit such an analogy further by representing the solution of 
the initial value problem explicitly in terms of the ultradiscretization (UD) of the 
Riemann theta function (z e M^): 

9(z) = fim elog ( expf — ^^^1^/2 + nz \ \ 

(1.1) ^ ^ V 

= - min|*nyln/2 + *nz}. 

Here A is the symmetric positive definite gx g integer matrix (|2.5(l appearing in the 
string center equation H4.1|l introduced in ^ . Likewise the Riemann theta function, 
0(z) enjoys the quasi-periodicity: 

(1.2) e(z + v) = *vA"i(z + v/2) + e(z) for any V G r = 

Let CL(n) = *nAn/2 -I- *nz be the quadratic form appearing in where L 

denotes the system size that enters A and z in our main formula H3.8|l . The ultra- 
discrete Riemann theta function Q{z) can be spotted in the following degeneration 

1 



scheme: 



exp(-CL(n)/e) 



n6Z9 



L ^ CO y 
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(1.3) 



min{cL(n)} = e(z) 



e{o,i}s 



UD \ 



^/ L ^ oo 



min |c(n)| 

ne{oa}9 



At the top there is the Riemann theta function, which degenerates into various 
objects. The UD procedure for getting 0(z) is the SE arrow from the top. 

Then in the hmit L oo, the minimum over n € shrinks down to that over n G 
{0,1}^, which reduces Ci(n) to its L-independent part c(n). Consequently, 8(z) 
tends to the bottom one in H1.3II , which we call the uUradiscrete tau function. The 
resulting expression H3.2|l for the infinite system gives the piecewise linear formula 
for the Kerov-Kirillov-Reshetikhin (KKR) bijection 8_ from rigged configurations to 
highest paths. One may go down the diagram H1.3|l via the other route. The thereby 
encountered function in the middle left is the sum of 2^ "trigonometric terms" that 
are characteristic in the tau functions of soliton solutions for the infinite system [Jj . 
In fact a procedure analogous to the SW arrow from the top has been described in 
p3.253 in jllj, where quasi-periodic soliton solutions tend to those in the infinite 
system. 

In our approach, the ultradiscrete Riemann theta function 0(z) arises most 
naturally by going from the bottom in (|1.3(l into the NE direction. The essential 
idea |l[)j is to embed a state p of the periodic box-ball system into an infinite system 
as p (g) p (8) p (g) • • • . It turns out that the ultradiscrete tau function for such periodic 
states is nothing but 0(z) up to irrelevant contributions. As an application we 
extend the problem to (C^)®^ and construct joint eigenvectors of the commuting 
time evolutions. The result may be viewed as an explicit formula of the Bethe 
vectors at gr = in terms of the ultradiscrete Riemann theta function. 

In section 121 we recall the periodic box-ball system and the inverse scattering 
algorithm that solves the initial value problem |ll)j . Section |21 contains our main 
theorem l3.3l Section^lgives the discussion on the connection with the Bethe ansatz 
at g = 0. 

We did not intend to make the paper completely self-contained. Exposition of 
the KKR bijection |Hj and Lemma rOl have been attributed to ^Hl- Rather, we have 
employed a casual description to clarify how the algorithmic solution to the initial 
value problem ^U] leads directly to the explicit formula (|3.8|l . We shall exclusively 
consider the case where the amplitudes of the solitons are all distinct, which greatly 
simplifies the presentation. The general case can be treated with the same idea. 

2. Periodic box-ball system and inverse scattering transform 

Let us quickly recall the periodic box-ball system without getting much into the 
crystal base theory. For a comprehensive treatment, see For a positive integer 
I, let Bi = {{xi,X2) & (Z>o)^ \ xi + X2 = 1} and set u/ = (^,0) £ Bi. The two 
elements (1, 0) and (0, 1) in Bi will be denoted by 1 and 2 for short. (Thus ui = 1.) 
In the following, the symbol ® meaning the tensor product of crystals can just be 
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{xi,X2) 
{Xi,X2) ® 2 1-^ 



understood as a product of sets. Define the map R : Bi ® Bi Bi ® Bi hy 

if (a;i,a;2) = (/,0) 
2 (g) + 1, a;2 — 1) otherwise, 

'2(g) (0,0 if (xi,X2) = (0,0 

1 ® {xi — 1, a;2 + 1) otherwise. 

i? is a bijection and called the combinatorial R. We write the relation R{u ®h) — 
b' (8) u' simply as w (g) & ~ h' ® u' , and similarly for any consequent relation of the 
form a®u®h®c:^a®h'®u'®c. 

A state of the periodic box-ball system is an array of 1 and 2, which is regarded 
as an element 6i g) • • • g) 6^ e Bf^ with L being the system size. Let the number 
of 2 S -Bi appearing in 6i ® • • • (g 6^ be Af. Without loss of generality we assume 
L > 2M (see ^Hl, section 3.3). Let V be the set of such states. Then the time 
evolution T; : P — > P is defined by 

(2.1) ui®p~p*®vi, vi® p ~Ti{p) ® vi. 

In the first relation, one applies the combinatorial R for L times to carry ui through 
p £ V to the right. This determines vi £ Bi and p* E V uniquely, {p* does not 
play an essential role.) Then the second relation using the so obtained vi specifies 
Ti{p), where the appearance of the same vi in the right hand side is a non-trivial 
claim (Unii section 2.2). vi is dependent on p as opposed to ui. 

The combinatorial R is the identity map on Bi® Bi, and therefore Ti is just the 
cyclic shift Ti (foi (g) • • • (g) &l) = &l (g 6i (g • • • (g 6l-i • The commutativity T;Tfc = TkTi 
holds for any k,l ([El, Theorem 2.2). 

Example 2.1. The time evolutions p,Ti{p), . . . ,T^{p) of the state p on the top 
line are listed downward for / = 2 and 3. The system size is i = 14. We omit the 
symbol (g. 

evolution under T2 evolution under Ta 

112 1112 2 2 1112 2 112 1112 2 2 1112 2 

2 2 12 11112 2 2 111 2 2 12 111112 2 2 11 

112 12 2 11112 2 2 1 112 12 2 2 111112 2 

2 112 112 2 11112 2 2 2 12 1112 2 2 1111 

2 2 2 12 1112 2 1111 112 12 2 11112 2 2 1 

112 2 12 2 1112 2 11 2 2 12 112 2 111112 

11112 12 2 2 1112 2 112 12 2 112 2 2 111 

2 2 1112 112 2 2 111 1112 112 2 1112 2 2 

112 2 112 1112 2 2 1 2 2 2 12 1112 2 1111 

2 1112 2 12 11112 2 1112 12 2 2 112 2 11 

Regarding 1 as an empty box and 2 as a ball, these patterns exhibit the nonlinear 
dynamics of balls. There are three solitons (wavepackets) with amplitudes 3, 2 and 
1 traveling to the right. 

Let us proceed to the direct and inverse scattering transforms. A state p — 
&i (g • • • (g &L is called highest if 

tt{l < i < fc I 6, = 1} > 11{1 < i < fc I 6; = 2} for all 1 < k < L. 

The state on the top line in example 12.11 is highest, whereas those on the second 
lines are not. Let V-^ be the subset of V consisting of the highest states. Any state 
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p E V can be expressed a.s p — Tf{p^) using some d G Z and a highest state p+ G 
7^+. For instance, the state T2{p) in example 12.11 is written as 22121111222111 = 
r^^(12111122211122). Given a state p, such d and p+ are not unique in general. 
Picking any one of them will be denoted by p i— > {d,p+). Consider the KKR 
bijection (f> from the highest state p+ to the rigged configuration ,8,: 



P+ 



«2 - 



(2.2) 

The partition (ig, . . . , 12, ii) is called the configuration and the integers < Ji < Pi 
are called the rigging. The combined data define a rigged configuration. Here pi is 
the vacancy number: 

(2.3) =L-2^min(i,j), 

where /i = {«i < 12 < • • • < *g}- Obviously, Pi^ > Pi2 > ■ ■ ■ > Pig holds, and it is 
known that ii + ■ ■ ■ + ig coincides with the number M of 6^ = 2 G -Bi contained 
in p+ = &i (g) • • • (g) 6l. Thus we have pi^ = L — 2M > by the assumption. See 
appendix A in jlOj for an exposition adapted to the present context. 

The configuration fi is actually independent of the non-uniqueness of the choice 
of and determined solely from p. The states are classified according to their 
configurations: 

where the disjoint union runs over all the partitions of M = 0, 1, . . . , [L/2]. V{ii) 
is the set of states whose configuration is Each subset 7^(/x) is invariant under 
any time evolution Tj, telling us that /i is a conserved quantity (jJOl Corollary 3.5). 
Physical meaning of /i is the soliton content, namely, the list of the amplitudes of 
the solitons involved in p. In particular g is the number of solitons. 

Unless otherwise stated, we shall consider those states whose configuration has 
the distinct parts as 

(2.4) ii = {ii<i2--- <ig}. 

Define the g y. g symmetric integer matrix A — {Aij)i,j,=^ and the lattice F by 

(2.5) = (5,jK + 2min(i, j), F = AZ^ c Z^. 

This matrix has arisen in the Bethe equation at q = H4.1() known as the string 
centre equation 9^. Under the condition L > 2M, A is positive definite. 

Let us proceed to the scattering data, i.e., the action-angle variables. The action 
variable is the set /i itself. The set of angle variables with prescribed p, is given by 
the quotient: 

(2.6) J = Jiii) - ZVF. 

The one to be assigned with the state p is found by the direct scattering map: 

$ : V{^i) ZxV+^ Jifi) 

^ ' ' p ^ {d,p+) ^ (J + dhi)/F, 
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where hi = (1, . . . , 1) e as defined in H2.8|l . J = {Ji)ii^^ G is specified by the 
KKR bijection as in H2.2|l . which we write as 4>{p+) = (/i, J) or simply (t>{p+) = J. 
Then J + dhi — {Ji + d)i^^. <i> is well-defined jn|. In particular, the non-uniqueness 
of the decomposition p {d,p+) is cancelled by taking mod F. For I e Z^, we 
denote its image in by the same symbol I. 

For I € we introduce the time evolution through 

(2.8) Ti(I)=I + h,, = (min(i,OW e ZS. 
Note that ihi = ^hi e F, therefore T/'(I) = Ie J. 

Theorem 2.2 (^Oli Theorems 3.11, 3.12). The map ^ is a bijection and the fol- 
lowing commutative diagram is valid: 

(2.9) T,| |t, 

Here Ti on the left and the right are given by \2.1]) and Y2.8\) . respectively. 

The composition o Tj o $ yields the algorithmic solution of the initial value 
problem by the inverse scattering method |S1 . The nonlinear dynamics Ti on 
V{ii) becomes the straight motion on J{n) with the velocity h;. In this sense 
J{^) is an ultradiscrete analogue of the Jacobi variety. Its cardinality is given by 
\J(pl)\ = detA = Lp^,p,,---pi^_, (Uni, (4.6), (4.13) and (4.21)). For I > ig, one 
has h; = hi^, hence Ti{p) = Ti^{p) by theorem 12.21 

In the limit L ^ oo, the quotient by F in (|2.6|l becomes void and the result 
provides the inverse scattering method for the box-ball system on the infinite lattice. 
The direct and the inverse scattering maps <I>^^ reduce to the KKR bijection (f>^^ 
itself. 

Example 2.3. For p ^ 22121111222111, let us derive 

(2.10) T^"°"{p) ^ 11112221112212, T^"°°{p) = 12211122111122 

based on the inverse scattering scheme (|2.9f) . (This p is T2{p) in example 12. II ') We 
have p = T^ip+) with the highest state p+ = 12111122211122. The image of the 
KKR bijection of (j){P+) ^'^d the direct scattering transform ^{p) are given by 



0(P4 







□ 1 






|3 







4 $(p) = 




2 
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Thus fi = {1,2,3}, ipi,P2,P3) = (8,4,2) and the matrix A (1^ reads 

fpi + 2 2 2 
A=l 2 P2 + 4 4 

\ 2 4 P3 + Q/ 

According to (|2.9|l and H2.8|l . the scattering data for the states T'2^3°"(p) are given 

by 
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rpWOO 



2003 



2006 
1002 



^1000 



($(p)) 



3003 



2006 



1002 



The angle variables appearing here are written as 

/1002\ /8\ /35\ /1002\ 

2006 ^ 4 + Ohi + ^ 161 , 2006 = 

\2003/ \1/ \161/ \3003/ 

The last terms involving A can be dropped by mod F, whereas the first terms in 
the right hand sides give rise to the rigged configurations and the corresponding 
highest states: 




11112221112212 ^ 











8 


4 



1 

11221111221221 







1 1 
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In view of +Ohi and +4hi , J^i^o" (p) and Tg^^o" (p) are obtained by taking the cyclic 
shifts r" and of these states respectively, in agreement with (|2.10|l . 

3. The explicit formula for the initial value problem 

First we present a piecewise linear formula for the KKR bijection. Let (/x, J) be 
a rigged configuration for a highest state in Bf^. To be concrete, we set 

r '((m, J)) = (1 - y(i), 2/(1)) v{L),y{L)) e V+, 

where y{k) S {0, 1} is the 'number of balls' in the k th box from the left. We 
parametrize the configuration ^ = {ii, . . . , ig} and the rigging J — {Ji^ , . . . , J^^) as 
in (|2.2I) . The following proposition 13 . 1 1 and lemma ^7]\ hold for the configurations 
such that ii < ■ ■ ■ < ig. 

Proposition 3.1. The image of the KKR bijection is given by 

(3.1) y(fc) =ro(A:)-To(fc-l)-Ti(fc) + ri(fc-l), 

(3.2) Tr(fc) = — min {> {Ji + ri — k)ni + > min{i, j)ninj} (r = 0, 1), 

ne{o.i}s — ' ^ — ' 

where n = (n^i, . . . ,nij. 

The proof will be given elsewhere for a more general case. Tr(k) E Z>o is the 
ultradiscrete tau function mentioned in section ^ We remark that there is no 
dependence on L in H3.2|l except in the upper bound pi (|2.3I) of the rigging Ji < pi. 
For k < 1 or k > L, (|3.1() gives y{k) = 0. As it turns out, after theorem 13.31 
proposition 13. II essentially provides the solution of the initial value problem of the 
box-ball system on the infinite lattice fc € Z. 

Lemma 3.2 ([13], Lemma C.l). Let q e Bf^ and r G Bf^ be the highest states 
associated with the rigged configurations (j^iq) — (A, I) and (f>{r) ~ (/i, J)- Then the 
rigged configuration of the highest state qi^r E j^®k+l (jji^q^^r) — (AU/i,IU J'), 
where 3' — {Jj)jefj. is given by 

J'j = J J +P]^ = - 2 ^ min(j, k). 

kex 



The shift pj here is nothing but the vacancy number in the rigged configuration 
(j>{q)- The notation {X(J fi,lU J') means the union regarding (A, I) and {n,J') as 
multi-sets of parts (rows in Young diagrams) assigned with rigging. For example, 



(A, I) 



(m,J') 



(AU/Lt,IU J') 







\ c 






a 




d 






b 





where, as usual, the ordering of the rigging d and b within a block of equal length 
rows does not matter. In what follows, we employ the convention of always arrang- 
ing the rigging to weakly increase upward within such blocks. 

Given a state p G V, take a highest state p+ e 7'+ and < d < L such 
that p = Ti{p+). Let (j){p+) — (a*, J) be the rigged configuration for p^, which we 
parametrize as fi = {ii, . . . ,ig} and J — {Ji-^, . . . , J^^). Here we assume ii < ■ ■ ■ < ig 
in accordance with the assumption (|2.4|l . We form a large highest state pf^ — 
p+ (g) • • • (E)p+ e Bf^^. By lemmaE21 its rigged configuration (^^, J^) (l){pf^) 
is given by 



3^ 



{il,l, ■ ■ ■ ,'il,N,i2,l, ■ ■ ■ ,i2,N, ■ ■ ■ ,ig,l, ■ ■ ■ , ig,N}, 

is, = + {a- (1 < a < N), 



where Pi = i — 2 ^^^^ min(i, j) is the vacancy number for p+ . We apply proposition 
IH.ll to (/i^, J^). From (|3.2|l the corresponding ultradiscrete tau function Tr{k) reads 

(3.3) 



™'^ivg'iX! X! ("^^a + ~ + X! X! min(i, j)ni,Qnj,^ 

ne{o,i}"9 I ije^ i<Q,/3<Ar 

where n = (nj^a, . . .^u^^^n, ■ ■ ■ I'^ig.i, ■ ■ ■,ni^,N)- Since Ji,i < Jj,2 < • ■ ■ < Jj,jv for 
each i (z fi^ the minimum here can be restricted to those n having the form 

— "-1,2 = • • • = rii^mi = 1, "■i,mi + l = n-i,mi+2 = • • • = Ui^jy = 

for some < < A^. Then the sums over a and /3 in (|3.3() can be taken, leading 
to 

E/ mi(mi — 1) \ ■r-^ . 

I mjJi H Pi + rriiri - niik ) + ^ mm[i,j)mimj 



(3.4) 



*m(j - I + rhoo - fchi) + i*mAm, 



where A = {Aij) is defined in (|2.5|l . We have set m = (mjig^, p = {Pi)iei-i ^nd 
used the vector notation J,hi,hoo around (|2.7|l and (j2.8(l . For instance hoo — h^^ 
and H2.3|l is rephrased as 



(3.5) 



p = Lhi - 2 ^ hj 



By taking iV to be even and shifting m to m+ hi, (|3.4|) is rewritten as *m(J— ^ + 
rhoo - (A: - ^)hi) + i*mAm + X, where X = f *hi(J- f +rhoo - {k-^)hi). 



This X can be put outside min, after which its dependence on r, k is cancelled in 
the difference (EH). Therefore we find that pf^ = (1 - y(l),2;(l)) ® ■ ■ ■ ® (I - 
y{NL),y{NL)) is given by H3.1|l with Tr{k) replaced by 

(3.6) Tr{k) = - min{*m(J - ^ + rh^ - {k - ^)hi) + i*mAm}, 

m Z Z Z 

where min is taken over those m — {mi)i^f^ E TP such that —N/2 < nii < N/2. 

From the relation^ = T{^(p+), the statep = (1 — a;(l), a;(l))(g)- • •0(1— a:(L), x{L)) 
is obtained from by picking up the length L segment corresponding to y{wL — 
d+ 1), . . . , y{{'w + 1)L — d) for any 1 < w < — 1. Thus in H3.6|) we replace k by 
k + wL — d with the choice w = to get Tr{k) = — mini„{cL(m)} with 

(3.7) CL(m) = *m(I - | - fchi + rh..^) + i*m^m. 

Here we have let I = J + dhi denote the angle variable <i>(p) for p. See (|2.7|l . 
The resulting formula for x{k) gives the state p corresponding to its action-angle 
variable (/x, I) as long as < d < L — 1, 1 < k < L and < Ji < Pi since we have 
started from the rigged configuration. These constraints are removed by taking the 
limit — > oo, where the minimum extends over m € Z^; therefore one has 

T,{k) = Q{l-^-kh,+rh^). 

By virtue of the quasi-periodicity of the ultradiscrete Riemann theta function H1.2|l . 
the difference 

x{k) = e(i - ^ - fchi) - e(i - ^ - (fc-i)hi) 

(3 8) 

- e(i - 1 - fchi + hoc) + e(i - 1 - (fc- i)hi + hoc) 

gains the invariance under k ^ k + L and I ^ I + v for any v G F = AZ^ . (Note 
that Lhi = ^hi £ F.) Namely, makes sense for k e Zl and I e J = Z^/r. 
To summarize, we have proved 

Theorem 3.3. For any state p G V of the periodic box-ball system, let (fj,, I) = *&(p) 
be the action-angle variable. Fix p = iPi)iefi by 51) and the matrix A by f2.5p . 
Then the state p is expressed as p ~ (1 — a;(l),a;(l)) (E) ■ • • (8) (1 — x{L), x{L)) with 
x{k) g {0, 1} given by j5'. Due to theorem \2.'/A this solves the initial value 
problem in that any time evolution T/'^ ■ ■ ■ T//* {p) is obtained by replacing I in 
with I + 7ih;i H h 7th;j (7.; e Z). 

The quadratic form (|3.7|l is decomposed as CL(m) = Lj^i^i rni{mi — l) /2+c{m.), 
where c(m) is independent of the system size L. In the limit L — > 00, the minimum 
is restricted to m e {0, 1}^ and O degenerates into the ultradiscrete tau function 
as in the scheme (|1.3|l . If I is chosen to be a rigged configuration, the formula (|3.8|) 
under such a reduction still describes the image of the KKR bijection although the 
function — minm6{o.i}s{c(m)} takes slightly different form from (|3.2() . The result 
provides the solution of the initial value problem of the box-ball system on the 
infinite lattice. 

In Figure 1, we plot the following function on the (fc,t) (space-time) plane: 



(3.9) u(fc,t) 



z?(T4(I) - E - fchi)^(r^(I) - E - (fc-l)hi + 
^(T^(I) - f - (fc-l)hi)^(T4,(I) - E _ fchi -Hh. 



where T^(T) = I + thoo by lES)) and i?(z) = Enezs exp(^-(*n>ln/2 + *nz)/ej is the 
Riemann theta function. In view of the scheme one has hnie^+o e log 0) = 
a;(fc). Thus u{k,t) gives a softening of the envelop of ultradiscrete solitons in the 
periodic box-ball system at e = under the time evolution Too- The selected 
parameters are 

For the periodic box-ball system described by H3.8() , this data corresponds to p = 
1122111111222222(g) 1®154^ which is a two soliton state with amplitudes 2 and 6. 
At t = 70, it becomes T^{p) = 1®*^"* (g 222222 (g l^^s ^ 22 l^^o. 
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J 



Figure 1. The envelope of the function u(k,t) (|3.9|l for 1 < fc < 
170, < t < 70. The top and right corners correspond to (fc, t) = 
(0,0), (170,0), respectively. It is periodic in the fc-direction. 

4. Discussion 

Theorem 13.31 enables one to construct the joint eigenvectors of Ti,T2,... in 
(C^)®^. The result may be regarded as an explicit formula for q = Q Bethe vectors 
in terms of the ultradiscrete Riemann theta function. We continue assuming that 
^ consists of distinct parts as in (|2.4|) . 

The Bethe equation for the periodic XXZ chain on (C^)®^ associated with 
Uq{sl2) becomes linear ai q = under the string hypothesis. The result is known 
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as the string centre equation [^: 

(4.1) Au=-| modZf, 

where u = (w^^ , . . . ,Ui^) with Ui being the centre of the length i string. We call u 
the Bethe root. In this normalization, the Bethe wave function is a rational function 
of exp(27rV— lui)i hence u lives in (M/Z)^. Thus there is one to one correspondence 
between the Bethe root u and the angle variable J G = V j AlP via the relation 

Au = J P. 

2 

The time evolution T/ of J 1)2. 8|l induces that of the Bethe roots, which is again a 
straight motion T/(u) = u + ^"^h; in (M/Z)f . 

At first sight, this appears contradictory, because Ti,T2, . . . are fusion transfer 
matrices at q = 0, which should leave the q — Bethe vectors invariant up to an 
overall scalar as well as the relevant Bethe roots. The answer to this puzzle is that 
the state p £ Bf^ that we are associating to u or J by <&(p) — (/i, J) is a monomial 
in (C^)®'^, which is not a Bethe vector at = in general. 

It is easy to remedy this. In fact, for each Bethe root u or equivalently J — 
Au+ ^ & J, one can construct a vector |J) S (C^)*^^ that possesses every aspect 
as a g = Bethe vector as follows: 

(4.2) |J) = ^ci,jp(I), 

cj = =p (-2,v/^'l(,4-'(J -R) + h.)y 

.(.) = (';,^y')«-«(';,l<^')eP(.,£(c»r. 

where x{k) e {0,1} is specified by We embed Bf^ into (C^)®^ naturally 

and extend Ti to the latter by C-linearity. The vector p{T) here is nothing but 
the state of the periodic box-ball system appearing in theorem 13.31 It follows that 
Ti{p{T)) = p{I + h(). Thus from J' + hi ^ J' , it is elementary to check 

THJ)=A,(J)|J>, 

Ai(J) =c-h„j =cxp (^2^yZTth,(u+^^)^ . 

The quantity A; (J) here exactly coincides with the q = Bethe eigenvalue given in 
equation (4.28) of jJUj. Note further that the transition relation 14.211 is inverted as 

where ci^j denotes the complex conjugate of ci j. It follows that the space of the 
q = Bethe vectors |J) coincides with the space of periodic box-ball states p for 
each prescribed soliton content /i, namely, 

C|J)= Cp. 
JeJ(M) pevifj.) 
Thus we conclude that the approach here bypasses the formidable task of com- 
puting the q ^ limit of the Bethe vectors in general, but leads to the joint 
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eigenvectors |J) of {Ti}. They form a basis of the space having the prescribed 
sohton content and possess the spectrum A; (J) anticipated from the Bethe ansatz 
at g' = 0. Moreover |J) is parametrized exphcitly in terms of the ultradiscrete 
Riemann theta function. 
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